Formation of a Direct Kolmogorov-like Cascade 
of Second Sound Waves in He II 
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Based on measurements of nonlinear second sound resonances in a high-quality resonator, we 
have observed a steady-state wave energy cascade in He II involving a flux of energy through the 
spectral range towards high frequencies. We show that the energy balance in the wave system is 
nonlocal in K-space and that the frequency scales of energy pumping and dissipation are widely 
separated. The wave amplitude distribution follows a power law over a wide range of frequencies. 
Numerical computations yield results in agreement with the experimental observations. We suggest 
that second sound cascades of this kind may be useful for model studies of acoustic turbulence. 
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Turbulence in the system of waves in a nonlinear non- 
dispersive medium is usually referred to as acoustic tur- 
bulence, or Burgers turbulence (BT) in honor of the 
Burgers equation which provides a useful description of 
the phenomenon. BT has been at the focus of numer- 
ous investigations during the last few decades because of 
its importance for basic nonlinear physics and in view 
of numerous applications in engineering and fundamen- 
tal science [l|, y, Well known examples of BT in- 
clude the turbulence of sound waves in oceanic waveg- 
uides [4], magnetic turbulence in interstellar gases [5], 
and shock waves in the solar wind and their coupling 
with the Earth's magnetosphere |6i]. 

In this Letter we propose a novel approach to the 
experimental study of BT: the use of nonlinear second 
sound standing waves in superfluid "^He (He II) within 
a high-quality resonator. Second sound is a slightly 
dissipative temperature/entropy wave that arises in su- 
perfluids and perfect crystals [7]. Second sound in 
He II is characterized by nonlinear properties that are 
rather strong compared to those of ordinary sound (pres- 
sure/ densitywaves) in conventional solids, liquids and 
eases [8,'9l,ll0llTT|. The velocity U2 of second sound de- 
pends on its amplitude and, to a first approximation, can 
be written as 



U2 = U2o{l + a5T), 



(1) 



where U20 is the wave velocity at negligibly small ampli- 
tude, ST is the wave amplitude, and 
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is the nonlinearity coefficient of second sound. C is the 
heat capacity per unit mass of liquid helium at constant 
pressure, and T is the temperature. 

There are huge advantages in the use of roton second 
sound waves in He II for model studies of nonlinear wave 



interactions. Within the experimentally convenient tem- 
perature range 2.17-1.5 K the nonlinearity coefficient a 
can be tuned, just by changing the bath temperature: 
a ^ — 00 near the normal-to-superfluid transition at 
Tx = 2.17 K; a passes through zero at — 1.88 K; 
and a - +2 K"^ at T - 1.5 K [8]. Thus one can 
study the dynamics of both nearly-linear and strongly- 
nonlinear waves with positive (like conventional sound) 
or negative nonlinearity while using exactly the same ex- 
perimental techniques. Such possibilities are unavailable 
in conventional experiments. The fact that the velocity 
of second sound iX20 ^ 20 m/s is more than an order of 
magnitude less than the velocity of conventional sound 
in gases and in condensed media allows us to increase the 
time resolution of the measurements. 

We use a thin-film heater to generate the second sound, 
and a thin-film thermometer (a low-thermal-inertial su- 
perconducting bolometer) as a detector. Use of a high- 
Q resonator enables us to create nonlinear standing sec- 
ond sound waves of high amplitude with only small heat 
production at the source. The wave amplitude ST can 
be changed from 0.05 mK up to a few mK, so that 
the Mach number can reach M = olST ^ 10~^ — 10~^, 
and the proper Reynolds number, defined [lH as Re = 
OLU2{){dST I dx) I ~ OiQST ^ can be changed from 1 to 
90 (here 7^; is the damping coefficient of a second sound 
wave of frequency cj, evaluated from the Q-factor of the 
resonator). It allows us to study, for the first time, the 
transition from the linear regime to turbulence in the sys- 
tem of second sound waves as the driving amplitude is 
increased. The experimental arrangements were similar 
to those used earlier [l^ . The resonator was formed by a 
cylindrical quartz tube of length L = 7 cm and inner di- 
ameter D = l.b cm. The film heater and bolometer were 
deposited on the surfaces of flat glass plates capping the 
ends of the tube. The heater was driven by an exter- 
nal sinusoidal voltage generator in the frequency range 
between 0.1 and 100 kHz: The frequency of the second 
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FIG. 2: Second sound power spectral amplitudes Af cal- 
culated numerically from Eq. (2) for four different driving 
force amplitudes: Fd = 0.01 (triangles); 0.05 (diamonds); 0.1 
(circles); and 0.3 (squares). The dashed line corresponds to 
Af oc as predicted by the theory [2, 18]. The arrow marks 
the boundary frequency fh of the inertial range for Fd — 0.3. 
The inset shows the calculated dependence on Fd of the stand- 
ing wave amplitude A: solid line - nonlinear waves, a < 
(T > Ta); dashed line - linear waves, a — ^ {T — To). 



FIG. 1: Power spectral amplitudes Af oi standing waves 
recorded at T = 2.079 K when driving at the 31st resonant 
frequency, fd = 3130 Hz. The ac heat flux density from the 
heater was (a) W 5.5 and (b) 22 mW/cm^ The dashed 
line in (a) is a guide to the eye; that in (b) corresponds to 
Af (X f~^'^ . The arrows indicate positions of the viscous cut- 
off frequency fb in each case. In the inset triangles show the 
observed dependence of the wave amplitude A on ac heat flux 
density W, measured for 31st resonance, and the dashed line 
corresponds to the Unear dependence A (x W valid at small 
W. 



sound (at twice the frequency of the voltage generator) 
was set equal to one of the longitudinal resonant frequen- 
cies. The measurements reported below were conducted 
at T = 2.079 K where a ^ -7.6 K-\ and at T = 1.775 
K where a ~ 1.8 K~^. The Q-factor of the resonator 
determined from the widths of longitudinal resonances 
at small heat flux densities was Q ~ 1000 for resonant 
numbers p < 10 and Q - 3000 for 30 < p < 100. The 
steady-state second sound distribution in the resonator 
was Fourier-analysed and its power spectrum was com- 
puted. 

Fig. 1 shows some typical experimental results at 
2.079 K. For a small ac heat flux density W < 1 
mW/cm^, corresponding to a standing wave amplitude 
less than 0.5 mK and to Reynolds numbers Re ~ 1, we 
observed a nearly linear regime of wave generation: the 
amplitude A of the standing wave was proportional to 
the heat flux density W (see Fig. 1 inset). Increase of 
the excitation above a few mW/cm^ (Re > 10) led to 
large deviations from A (xW and to visible deformation 



of the initially sinusoidal standing wave, accounted for 
by the formation of multiple harmonics in its spectrum, 
as shown in Fig. 1 and the inset. At W ~ 80 mW/cm^ 
we observed a marked reduction of the Q-factor to ~ 500 
for longitudinal resonances with p ^ 30. This may be 
associated with e.g. broadening of the resonance curve 
due to vortex creation in the bulk He II. In what fol- 
lows, we discuss only the measurements made at < 80 
mW/cm^. 

It is clearly evident from Fig. [1] that the main spec- 
tral peak lies at the driving frequency fd, and that high- 
frequency peaks appear at its harmonics, fn = fd x ^ 
with n = 2,3,... It can be seen that a cascade of waves 
is formed over the frequency range up to 15 kHz in 
Fig. Ufa), and up to 80 kHz in Fig. [IJb), i.e. up to 
frequencies 5 and 25 times higher than the driving fre- 
quency. As shown in Fig. [IJb), the dependence of peak 
height on frequency may be described by a power-like law 
Af = const X f~^ at frequencies lower than some cut- 
off frequency f^. For sufficiently high heat flux densties 
W > 10 mW/cm^ the scaling index tends to s ~ 1.5 ±0.3. 

In order to understand how the wave cascade is formed, 
we have undertaken a numerical study of the dynamics 
of nonlinear second sound waves within the high-Q res- 
onator. Full details will be given elsewhere, but in essence 
it involves direct integration of the two-fluid thermohy- 
drodynamical equations [7], expanded up to quadratic 
terms in the wave amplitude. We use a representation of 
the second sound waves in terms of Hamiltonian variables 
[H, [l4|, permitting us to improve the precision of the 



calculations. The numerical technique is similar to that 
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used in earlier studies [H of developed two-dimensional 
acoustic turbulence. In the present model, however, wave 
damping was taken into account at all frequencies, a fea- 
ture that is of key importance [15] for a correct descrip- 
tion of the formation of a cascade of nonlinear sound 
waves with increasing driving force amplitude. The nu- 
merical results obtained are suitable for direct compari- 
son with our experimental observations shown in Fig. 1. 

The equations of motion governing non-local energy 
balance in the system read 



^ — ^ ^ ^n,n\,n2Q^n-i^n'2^Tt—n\—n2 

ni ,77,2 



i) - ilnK + Fd, 



(2) 



where hn{t) = {l/2){B-^Sn + iBnPn) is the time- 
dependent canonical amplitude of second sound at the 
n-th resonant frequency fn] Sn and /3n are the space 
Fourier components ~ cos{27rfnx/u2oL) of the entropy 
and of the potential of the normal-superfluid relative ve- 
locity; Bn = (/nC/T)^/2; Vn,ni,n2 ^ a(nnin2)^/^ de- 
scribes the three-wave interaction; 7^ = un'^ models the 
viscous damping of second sound; (x W is the am- 
plitude of the force driving the n-th resonant mode. We 
neglect the coupling between temperature and pressure 
oscillations in view of the smallness of the thermal ex- 
pansion of He II in the relevant temperature range. The 
wave spectrum was calculated as cx Bn {bn + 6*). We 
used a periodic driving force of frequency equal to a res- 
onant frequency of the resonator corresponding to the 
conditions of the measurements, taken as fd = 3000 in 
dimensionless units. The effective viscosity coefficient u 
was fitted to the measured value of the quality factor 
Q ~ 3 X 10^ of our resonator, to facilitate comparison of 
the model results with those from the experiments. 

Fig.[2]shows the evolution of the steady-state spectrum 
with increasing driving force amplitude Fd^ calculated for 
a < {T > Ta). Points on the plot correspond to the 
amplitudes of the peaks in the spectrum. It is seen that, 
at small driving amplitude Fd ~ 0.01 (triangles), viscous 
damping prevails at all frequencies and a turbulent cas- 
cade is not formed: the amplitude of the second harmonic 
is an order less that the amplitude of the main harmonic. 
In this regime the wave shape is close to linear. At inter- 
mediate driving amplitude ~ 0.05 (diamonds) nonlin- 
earity starts to play a role at frequencies of the order of 
driving frequency, and a few harmonics are generated (cf. 
Fig. 1(a)). At high driving amplitudes Fd > 0.1 (circles 
and squares) a well-developed cascade of second sound 
waves is formed up to frequencies 30 times higher than 
the driving frequency, i.e. the behaviour predicted by the 
computations is very similar to that observed in our ex- 
periments (Fig. 1(b)). 

Formation of the spectra observed in the experiment 
and numerics is evidently attributable to the cascade 
transfer of wave energy through the frequency scales due 



to nonlinearity, thus establishing an energy flux in K- 
space directed from the driving frequency towards the 
high frequency domain^ Following the basic ideas formu- 
lated in Refs. P, [H, 



17 



18[, such a highly excited state 
of a system with numerous degrees of freedom is referred 
to as turbulent. So we may infer that at relatively high 
driving amplitudes we observe acoustic turbulence formed 
in the system of second sound waves within the inertial 
(non-dissipative) range of frequencies. Formation of the 
observed direct cascade is qualitatively similar to cre- 
ation of the Kolmogorov distribution of fluid velocities 
over frequencies in the bulk of classical inconipressible 
fluids [7], or in a system of surface waves [H, 0. But 
we found both in experiment and numerics that the am- 
plitudes of high frequency harmonics rise gradually with 
the increase of the driving amplitude. 

In some range of frequencies the calculated spectrum 
obeys a power law with 5 ~ 1 (squares in Fi g. 2 ) that 
is in qualitative agreement with the theory fj. Iisj. The 
difference in slope of the calculated spectrum from the 
value s = 1.5 ± 0.3 observed in the experiment may 
be attributable to our resonator having an insufficiently 
high quality factor, as well as to the influence of other 
mechanisms of wave energy relaxation not taken into ac- 
count in our computations. The precise identification of 
the scaling index characterising a turbulent distribution 
is obviously important for connecting BT theory to real 
physical experiments. To clarify the origin of the differ- 
ence between observations and computations, we plan to 
make further measurements and to undertake a system- 
atic study in the vicinity of Ta, where the nonlinearity of 
second sound rises much faster than its damping [11], and 
the proper Mach and Reynolds numbers should therefore 
be substantially increased. 

From a comparison of the measurements and computa- 
tions (Figs. [T]and[2l insets) we can conclude also that the 
deviations from the linear dependence A (x W observed 
at moderate driving amplitudes is caused by nonlinear 
outflow of the energy to higher frequency scales, eventu- 
ally reaching the dissipative frequency range, rather than 
by the scattering of second sound from quantized vortices 
in He II, as considered earlier [III, [2^ . 

The effect of the high frequency cutoff of the wave 
spectrum manifests itself as an abrupt decrease in the 
amplitudes of the harmonics at W < 10 mW/cm^ (com- 
pare Fig.[T](a) and Fig. [2] (triangles)), and as a change of 
slope on double- log scales, at higher W (cf. Fig. [TJb) and 
Fig. [2] (circles and squares)), at some frequency f = fb- 
f ^ fh the nonlinear mechanism for nearly non- 
dissipative transfer of the wave energy changes to vis- 
cous damping of the waves (cf. observations of the high- 
frequency edge of the inertial range of frequencies of cap- 
illary turbulence on the surface of liquid hydrogen [l5[ ) . 
It causes a faster reduction of sound amplitude at fre- 
quencies / > as observed. 

The dependence of the boundary frequency on the 
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in qualitative agreement with the results of our experi- 
ments. 
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FIG. 3: Dependence of the viscous cutoff frequency fh on the 
amplitude A of the standing wave for different temperatures, 
for resonance numbers p = 31 and 32. Dashed lines: the 
results of our computations based on the theory [18]; data 
points: experimental measurements. 



standing wave amplitude A is shown in FigO It is clearly 
evident that the inertial range is extended towards higher 
frequencies when the driving amplitude is increased. At 
sufficiently large amplitudes A, the boundary frequency 
fh = const (T, fd) X A, which agrees with our calculations 
based on [18]. Such a dependence was observed when 
driving at resonant frequencies with odd resonant num- 
bers p (the full symbols in Fig. [3]). 

When driving at even p (open symbols), however, the 
boundary frequency is noticeably lower than that mea- 
sured for the nearest odd resonant number, for ly > 10 
mW/cm^. This reduction may be connected with a 
change in the mechanism of energy relaxation in the wave 
system caused by the generation of subharmonics with 
frequencies lower than fd- We recently detected the ap- 
pearance of subharmonics at high W when driving at 
even p, a phenomenon that promises to be of huge inter- 
est for nonlinear sound wave dynamics [12|, |2l| . We plan 
to study it in detail in the near future. 

In conclusion, we have demonstrated that the system 
of roton second sound waves in a high-Q resonator filled 
with He II can be used as an effective tool for the de- 
tailed modelling and investigation of acoustic turbulence. 
We observed for the first time a smooth crossover in the 
system of second sound waves from nearly-linear regime 
at low driving amplitudes (at Re ~ 1) to a nonlinear 
regime at moderate driving amplitudes, and, further, 
to developed turbulence at high driving amplitudes (a 
Kolmogorov-like cascade [l6| at Re ~ 10^), rather than 
a sharp transition to turbulence at Reynolds numbers 
higher than some critical value as occurs in incompress- 
ible fluids [7]. In the high frequency domain a cutoff of 
the cascade is observed, caused by a change in the mech- 
anism of energy transfer, from nonlinear wave transfor- 
mation to viscous damping. Numerical calculations are 
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